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Abstract 

The time evolution is studied for the Landau level problem with a general time dependent electric 
field E(t) in a plane perpendicular to the magnetic field. A general and explicit factorization of the time 
evolution operator is obtained with each factor having a clear physical interpretation. The factorization 
consists of a geometric factor (path-ordered magnetic translation), a dynamical factor generated by the 
usual time-independent Landau Hamiltonian, and a nonadiabatic factor that determines the transition 
probabilities among the Landau levels. Since the path-ordered magnetic translation and the nonadia- 
batic factor are, up to completely determined numerical phase factors, just ordinary exponentials whose 
exponents are explicitly expressible in terms of the canonical variables, all of the factors in the factor- 
ization are explicitly constructed. The numerical phase factors are quantum mechanical in nature and 
could be of significance in interference experiments. The factorization is unique from the point of view 
of the quantum adiabatic theorem and provides a demonstration of how the quantum adiabatic theorem 
(incorporating the Berry phase phenomenon) is realized when infinitely degenerate energy levels are in- 
volved. Since the factorization separates the effect caused by the electric field into a geometric factor 
and a nonadiabatic factor, it makes possible to calculate the nonadiabatic transition probabilities near 
the adiabatic limit. A formula for matrix elements that determines the mixing of the Landau levels for 
a general, non-adiabatic evolution is also provided by the factorization. 



1 Introduction 



The magnetic translation concept is much discussed in the literature [TJ [3J [31 IH [HI H] • The basic example is 
the quantum mechanics of a charged particle moving in a two-dimensional plane perpendicular to a magnetic 
field, i.e., the Landau level problem. The Hamiltonian does not commute with the usual translation operator 
generated by the canonical momentum. The translation symmetry of the physical situation is realized 
through magnetic translation. However, in this simple context, the solution of the problem (energy levels, 
etc.) does not seem to rely on this concept in any essential way. The magnetic translation concept has 
been used EJ [TO] [IT] in condensed matter physics including Bloch electrons in a magnetic field where 
there is a lattice potential present. Yet we believe the essential role of magnetic translation in exploring 
time-dependent problems have yet to be fully explored. 

In this paper, we study the quantum mechanical problem of a charged particle moving in a two- 
dimensional plane subjected to a uniform magnetic field perpendicular to the plane and a spatially uniform 
but time-dependent in-plane electric field E(t). The Hamiltonian for such a system is the sum of the usual 
Landau level Hamiltonian H and the potential energy of the charge in the electric field: 



where A(x) satisfies V x A(x) = Be^. The situation contains the usual Landau level problem as a special 
case, where Ei(t) = E2(t) = 0. It also contains the special case for which the uniform electric field is constant 
in time, and the case [12] for which E(i) depends on time but is along a fixed direction. These cases are 
all known in the literature. Our purpose here is to study the general case with E(t) = Ei(t)ei + E2(t)e2 
being a general time-dependent electric field. This general situation cannot be reduced to the sum of two 
special cases and it will be shown that it has new features of its own. In particular, we will show that a 
path-ordered magnetic translation, which is an element of the magnetic translation group, plays an essential 
role in describing the evolution of the system. This path-ordered magnetic translation contains a numerical 
phase factor that is nontrivial only when the electric field changes direction with time. Our main result is 
a factorization of the time-evolution operator into three factors, each having a clear physical meaning. The 
method and result seem to be quite natural that it is possible that they may find applications in studying 
more complicated situations such as the Landau problem on a cylinder [13j or a torus, or when there is a 
lattice potential present. 

First of all, from the physical picture that underlies the Hall effect, the circular orbit of a charged particle 
in mutually perpendicular electric and magnetic fields exhibits motion in the direction of E x B. It is easy to 
conjecture that, in the case where E(t) is a general time-dependent electric field, the instantaneous velocity 
of the global motion of the circular orbit is cE(t)/B, in the direction of E x B. Therefore, the position of 
the center of the circular orbit at time t is described by a parameter 



where we have taken the initial position R(0) to be at the origin of the parameter space: R(0) = 0. 

In this paper, we aim to find a general factorization of the time-evolution operator U(t, 0) corresponding 
to H in terms of three factors: a geometrical factor describing the physical displacement of a quantum wave 
by the amount of R(i), a dynamical factor generated by the usual time-independent Landau Hamiltonian 
Hg, and another factor that describes the mixing of the Landau levels of Hq. The factorization is general in 
the sense that it is valid for a general time variation of the function R(i) and therefore of E(i). 

We comment here that previous work on the special cases, such as the one in [12] has not relied on the 
magnetic translation concept. This is due to the fact that if the electric field is along a fixed direction, say 
in the x\ direction, one chooses the Landau gauge where A(x) = (0,Bx,0) and furthermore specify the 
constant of motion p2 to be and the problem is then simplified. Such a simplification also implies that the 
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— [p - -A(x)] 2 - qEttyxi - qE 2 (t)x2, 
2m c 



(1) 




(2) 



2 



quantum states under consideration have p 2 = 0. For such a restriction, the geometric operator found in 
the present paper is equal to the identity. So the only effects caused by the electric field are nonadiabatic 
transitions. No such simplification exists if the electric field changes direction with time. In order to treat 
the general case, the magnetic translation concept is necessary. In fact, even in the special case, if one does 
not make the special restrictions, there should be a magnetic translation along the fixed direction of X2- We 
also note that upon choosing a specific gauge, the Hamiltonian considered here can be cast in a quadratic 
form. For such type of time-dependent Hamiltonians, there is a method that uses the corresponding classical 
solutions to construct the propagator of the quantum problem [IH [15] . While such an approach could be 
useful as a general theory, it does not seem to point to a factorization of the time evolution operator which, 
in specific contexts, can make the time evolution transparent. 



2 Factorization of the time-evolution operator 

In this section, we first give a discussion on the properties of the variables (7ri,7r 2 ) and (101,102) which 
can be expressed in terms of the canonical variables x and p. Then we use these variables to construct 
a factorization of the time evolution operator. Although the commutation relations among these variables 
hold in any picture, we assume throughout this paper that when these variables appear in the exponents, 
they are Schrodingcr picture variables. 

In the usual Landau problem Hq, the kinematical momentum is 7r M = p^ — |A M (x), /i = 1,2, where A^x) 
is the vector potential in arbitrary gauge, though for simplicity we assume that ^(x) does not depend on 
time. Define 

qB qB qB qB 
W 1 = 7Tl x 2 = C 2 , W2 = TT 2 -\ £1 = — c 1 . (3) 

c c c c 

In classical mechanics, (01,02) is the center of the circular motion of the charged particle in the magnetic 
field. In quantum mechanics, because of the canonical commutation relations, c\ and c 2 do not commute. 
We have the following commutation relations 

[tti,^] = ihqB/c, [wi,w 2 ] = —ihqB/c, [n^, w u ] = 0. (4) 

To realize a translation x — > x + R(i), where R(i) traverses through a path Cr in parameter space, one may 
use the ordinary translation operator exp^—ip^R^ty/h). (Summation over repeated indices is assumed.) 
However, the ordinary translation does not commute with n^. The magnetic translation operator, defined 
as 

M(R) = exp(-iw M R M /h), (5) 

is a generalization of the ordinary translation operator when a magnetic field is present: It physically 
translates a quantum wave because it commutes with the kinematical momentum, and in the simplest 
Landau system Hq, it preserves the energy. Note that unlike ordinary translation, a distinction between 
exp^—iw^R^/h) and the path-ordered magnetic translation 

M(C R ) = Pcxpt-i/rV f dR^) = Pexpi-ih^wM (6) 

has to be made, for a general path Cr that is not in a straight line, because w\ and u>2 do not commute. 
(We assume i? M (0) = 0.) Because of the simple commutation relation [101,102] = —ihqB/c, their difference 
is a numerical phase factor, i.e. 

M(Cr) = e i ^ c ^exp(-iw fl R fl /h), (?) 

where /3(Cr) is a real number determined by the path Cr traversed by R(t). In particular, for a closed 
path, /3(Cr) is equal to — Ig, where is the magnetic flux enclosed by the loop C. This follows from the 
definition of a path-ordered exponential, the commutation relation [101,102] = —ihqB/c and the formula 
e A e B = e A+B e^ A,B ^ for A and B commuting with [A, B]. For an open path R(i), the flux is through the 
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area enclosed by the path and the straight line pointing from the end point R(t) to the initial point R(0) = 0. 
In general we have 

Or, if we denote R(t) = Ri(t) + iR 2 (t), then 

/3(C R ) = -^^/ (R*dR - RdR*). (9) 
hc 4 Jen 

Another path-ordered exponential that is relevant to our purpose is 

J(C U ) = Pexp(wu/h-^u*/h), (tt = TTi + in 2 ) (10) 

where the path C u is the one traversed by u(t) in a complex u-plane. Similar to the path-ordered magnetic 
translation, this path-ordered exponential can be evaluated by using the formula [n,^] = 2hqB/c. We have 

J(C U ) = e^ (c »>exp -TiV/ft), (11) 

where 

7 (C„) = 4^ / (u*du - udu*) = -^4S(C U ). (12) 

S(C U ) is the area enclosed by the path traversed by u(t) in the complex w-plane and the straight line 
connecting the end and initial points of the path. 

The basic idea that leads to the factorization of U(t, 0) is to switch from the Schrddinger picture to the 
Heisenberg picture first. Furthermore, we observe that the Heisenberg equations of motion decouple from 
each other for the variables 7r M (i) and w^(t), unlike the equations of motion for the canonical variables. 
It is the behavior of n^t) and w^(t) instead of that of and x^t) that leads to the factorization of 

U(t,0); namely we find an operator 0(t, 0) in a factorized form that recovers the evolution of 7r M (t) and 
w^t) through 7r M (t) = O^(t,0)TTfj,O(t,0) and w^t) = 0^(t, 0)w IJ ,O(t, 0). We then verify that this operator is 
in fact U (t, 0) by a computation that shows it satisfies the Schrddinger equation with the initial condition 
U(0, 0) = I. These steps are presented in detail in Appendix A. The exact and general factorization of the 
time evolution operator U (t, 0) is then determined to be 

U(t,0) = M(Cn)D(t)J(C u ), (13) 

where D(t) = cxp(— iHot/K) is the evolution generated by the usual Landau Hamiltonian; M(Cr) de- 
scribes path-ordered magnetic translation of a wave corresponding to the path Cr traversed by R(i) = 
(E 2 (s)ei — Ei(s)e 2 )ds. These two operators commute with each other, so no energy is gained or lost 
by the action on the wave-function of a magnetic translation. The complex parameter u that determines the 
operator J(C U ) is given by 

u(t) = i -j\- i "°j-R*{s)ds = Zlj\-^E*{s)ds, (14) 

where R*{a) = R^s) - iR 2 {a), E*(s) = E^s) - iE 2 {s), and uj = qB/(mc). 

We see that the operator exp (wu/h — ir^u* /tij in J(C U ) mixes different energy levels of H a . This is 
due to the fact that tt 1 ' and tt have the meaning of being proportional to the creation and annihilation 
operators: the Hamiltonian H can be written as H = Hco(a^a + 1/2), where = {hk)" 1 ^ ,a = (fifc) _1 7r, 
with (fifc) -1 = y/ c/(2qBh). From this we see that the operator J(C U ) has the meaning that, when acting on 
any of the ground states of Ho, it generates a coherent state associated with the minimization of Atti ■ An 2 
rather than the minimization of AxiApi, or of Ax 2 Ap 2 . The coherent state nature of J(C u )\^o), where |*o) 
is a ground state of Hq, is then preserved under the action of M(Cn)D(t) in the time-evolution operator. 
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3 Gauge invariance and explicitness of the factorization 



Apart from the numerical phase factors e 1 ^^' and e i ^^ Cu \ the operators M(Cr) and J(C U ) are generated 
by (wi,W2) = (— ^7-02, ^§-Ci) and (tti,^) respectively. In the system H , the physical meaning of these 
generators are given by the "center of the circular orbit" and the kinematical momentum which are gauge 
independent. However, one must be careful in claiming the factorization gauge invariant in the most general 
way since the derivation of the factorization relies on the condition that 7r M and in the exponents are 
time-independent Schrodinger variables and as such, is assumed to depend on x only. For example, in 
the symmetric gauge where A(x) = -j e 3 x x > we have m — pi + ^£2, "2 = P2 ~ ^f^i; w i = Pi ~ ^B x 2 and 
w-2 = Pi + %-X\. Therefore, the three factors in the factorization are all explicit functions of the canonical 
variables once a specific gauge is chosen. 

Consider the example of a rotating electric field. We have R(t) — R^(e~ tut — 1) or E(t) = (iB/c)R(t) = 
(vB/c)R e~ iut = E e~ tut . By the formulas (9), (12) and (14), we have /3(C R ) = {qB / hc)\R Q 2 {vt - sini/i), 
u(t) = {-R v/2){e- lu>t+lut -l)/(-iu + iu), and j(C u ) = (qB/hc){R Q 2 /2)( : ^) 2 [(uj - v)t - sin(w - v)t]. For 
the case of v — > u, we have u(t) = —R i/t/2, and j(C u ) — 0. The expressions for A/(Cr) and J(C U ) arc 
therefore obtained according to (7) and (11). In this particular example, the result implies that resonance 
happens at v = u) where u(t) increases linearly with time, thereby causing rapid transitions among the 
Landau levels. This example will be further commented on at the end of the paper. 

4 Physical implications of the factorization 

The general factorization of U(t, 0) demonstrates that a path-ordered magnetic translation is a natural 
concept associated with a charged particle in a time-dependent electric field and a uniform perpendicular 
magnetic field. The same can be said of the path ordered exponential J(C U ) that describes the mixing of 
the Landau levels. The factorization also has a natural connection with the quantum adiabatic theorem as 
discussed in the next section. The numerical phase factors e n<,Cu ^ in J(C U ) and e l ^ c ' R - ) in M(Cr,), which 
distinguish the path-ordered exponentials from the direct exponentials, represent pure quantum effects that 
have no classical origin and they could be of consequences in interference experiments. It is obvious that 
the factor e l ^ c ' R ) is nontrivial only when the electric field changes direction with time. This numerical 
factor contains the adiabatic Berry phase in the usual sense as discussed in the next section. One may want 
to draw an analogy with the Aharonov-Bohm phase even though the magnetic field interacts directly with 
the particle here. However, the phase /3(Cr) = — ^ that results from a path-ordered magnetic translation 
around a closed path Cr has an opposite sign from the Aharonov-Bohm phase. The derivation of the latter 
from the point of view of the geometric phase is given in Berry's paper [16] . 

5 Relation to the quantum adiabatic theorem 

The factorization can be viewed from a different though equivalent perspective. The Hamiltonian H can be 
recast in the form 




(15) 



through the gauge transformation 



*i(x,t) = expH|-x(x,R)]*(x,t), 



(16) 



A L (x, R) = A(x) - Vx(x, R) = A(x) - Be 3 x R(t) , 



(17) 



where 



X(x,t) 



BR 2 (t)xx +BR 1 {t)x 2 . 



(18) 
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Because ^(x, 0) = ^(x, 0), the time-evolution operators E/r,(i, 0) and U(t, 0) are related by 

U L (t, 0) = ex P H^ x (x, H)]U(t, 0). (19) 

The Landau level Hamiltonian unlike the gauge equivalent H(t), has energy eigenvalues E n — frw(n + 

1/2), where u) — qB/(mc), that depend on B only and is independent of the time variation of HL,{t). 
Therefore, it is in the gauge of HL(t) that D(t) carries the dynamical phase factor of adiabatic evolutions 
of eigenstates of the Hamiltonian. This should relate to the quantum adiabatic theorem, where G(Cr) = 
exp[— «^x( x j R)]-^(Cr) is a geometrical operator completely determined by the path of R(i) which brings 
an initial eigenstate of Hl(0) to an instantaneous eigenstate of H^lt), and J(C U ) describes nonadiabatic 
transitions. This is an example of the quantum adiabatic theorem |17j where all of the three factors of 
the time-evolution operator are explicitly constructed and where the energy levels of the instantaneous 
Hamiltonian are infinitely degenerate. 

Recall that the usual quantum adiabatic theorem (and also generalizations 18 )incorporating the Berry 
phase phenomenon is essentially a factorization of the time-evolution operator into three pieces: a geometric 
factor that embodies a Berry phase, a usual dynamical factor, and a nonadiabatic factor that approaches the 
identity operator in the adiabatic limit [H]. In our case, once a specific gauge is chosen, the factors G(Cr) = 
exp[— i^x( x J R)]-^(Gr) and J(C U ) are exponentials of explicit functions of the canonical variables, it thus 
provides an explicit example of the quantum adiabatic theorem involving infinitely degenerate energy levels. 

In a previous work [S], we obtained a factorization of the time evolution operator for a charged particle 
in a slowly rotating magnetic field with a strong confining potential confining the particle to be in the plane 
that is perpendicular to the instantaneous magnetic field. There, the factorization is valid in the adiabatic 
limit only; i.e., no information about the nonadiabatic factor was obtained due to the complexity of the 
problem. The method adopted in the present paper may be applied to study the nonadiabatic factor in that 
situation. 

6 Adiabatic perturbations 

It is clear that because of the existence of the factor M(Gr) in the time-evolution operator U(t, 0), which 
gives rise to the geometric phase phenomenon, one cannot choose a fixed basis of eigenfuctions of Ho, and 
perform a standard textbook version perturbative calculation. This is true even if the electric field is small, 
because it is the accumulative effects, i.e., R in M(Gr), that determines M(Cr) which is not close to identity 
even if E is small. 

The factorization of U(t, 0) allows us to do this perturbative calculation precisely because it identifies a 
set of parameter-dependent bases with the help of M(Cr), then nonadiabative transition probabilities are 
completely determined by J(C U ). 

The expression for J(C U ) therefore allows the explicit calculation of non-adiabatic transition probabilities. 
Take an initial eigenstate |^(n, R(0))) = |W(n, 0)) of the Hamiltonian Hl(0) with eigenvalue E n — huj(n + 
1/2). The factorization of Ul says that the time-evolution of |^(n, 0)) can be seen as the action of G(Gr)Z?(£) 
on top of J(C u )\Q(n, 0)). Since G(C-R,)D(t) does not cause transitions, the transition is caused by the action 
of J(C U ) on |$(n, 0)} only. In the adiabatic limit where |R(£)| ~ u(t) is small, one can expand J(C u )\Q(n, 0)) 
into 

J(G„)|$(n,0)) = e i7(c " ) (l-|-(7ru/ft-7r t u7^)|$(n,0))) + O(u 2 ), 
= e^ c "}(l + (uka - ku*ai)\$(n, 0))) + 0{u 2 ). 

The transition probabilities are 

|($(n-l,0)|J(C„)|$(7 1 ,0))| 2 = |($(7i-l,0)| U M$(n,0))| 2 + O( U 3 ), 

= nk 2 \u(t)\ 2 + 0(u 3 ), 
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K*(n + 1, 0)| J(C u )\t>(n, 0))| 2 = |($(n + 1, 0)|aV/(ftfc)|$(n, 0))| 2 + 0(u 3 ), 

= (n+l)(^)Kt)| 2 + 0( M 3 ). 

All other transition probabilities are zero to the order of 0(u 2 ). So we see in this case that non-adiabatic 
transition probabilities are dependent on the energy levels, and increase with n. This is a result that shows 
explicitly how the quantum adiabatic theorem can be realized when an infinitely degenerate energy level is 
involved. As usual, the role of the oscillating e" tus in the expression for u(t) is to make the effect of E*(s) not 
to accumulate during an adiabatic process with t £ [0, T], where T is a parameter that R and E may depend 
on through R(^) and E(^). If Eq/B is used as the dimensionless small parameter, say the electric field has 
a constant magnitude and that E changes direction slowly compared with u, then T can be chosen to be 
{B/E)u)~ x . When T > one estimates that \u(t)\ <~ cE/Buj, for t £ [ 0,T]. The transition probabilities 
are then bounded by n(cj~ 1 c/l B ) 2 • (E 2 /B 2 ) for t £ [0,T], where l B = \Jhc/{qB) is the magnetic length. 
For an electron in a 15 T magnetic field, if E = 1000 volts/m, we have E/B = 1000/(15 x 3 x 10 8 ), then 
n(u}- 1 c/l B ) 2 ■ (E 2 /B 2 ) = n ■ 1.45 x 10" 5 , for the duration of T = 1.71 x 10" 3 s. Then for a ground state for 
example, excitations can be expected to occur in about 100 s. 



7 Mixing of the Landau levels beyond the adiabatic limit 

One may also study the mixing of the Landau levels for a general nonadiabatic time evolution. Since 
M(Cn)D(t) in the time evolution operator does not cause transitions among Landau levels of Hq, the 
transitions are caused by the action of J{C U ) on |4>(n,0)) only. In the most general case, the expression of 
|($(m, 0)| J(C u (t))\$(n, 0))| 2 needs to be evaluated in order to determine the transition probability from an 
initial state (at t = 0) that is at the n-th Landau energy level to an m-th energy level at time t. We have 

J(C U ) = e^ c ^exp(mi/h-^u*/h), 
= e i ^ c ^exv(uka-u*k(J). 

Using the formula e A+B = e A e B e~^ A - B \ and the commutation relation [a, a*] = 1, we have 

J{C U ) = e *7(C„) e -i|«fc| 2 e -«*fcat e uka _ 

Therefore, we have the following general expression for the matrix elements: 

(m\J{C u (t))\n) = e^^e-^ 2 {e- uka \m))^ (e uka \n)), (20) 

which represents mixing of the Landau energy levels for the most general type of the electric field which is 
not necessarily small and which does not have to change slowly. 

In particular, this formula implies that for the ground state we have the following matrix element: 

(0\J(C u {t))\0) =e i ^ c «>e-5l« fe l 2 , 

which implies that if we start with a ground state of Ho at time 0, the probability that it remains to be 
in a ground state is (0| J(C u (t))\0)\ 2 = e~' ufe ' = e~l"l "kh". This implies, for the rotating electric field 
example as we discussed earlier where E — Eoe~ lult and u{t) = —cEot/(2B) at the resonance of v = co, the 
probability for the state to remain in a ground state of Hq (not necessarily the initial ground state because 
of the existence of the magnetic translation) is exactly 

P r (0 -> 0) = cx V [-2{E /B) 2 {c 2 t 2 /l B 2 )l (21) 

where l B = y/fre/qB is the magnetic length. (Note that for the electron, q = — e, the resonant electric field's 
angular velocity is then along the positive z direction, as expected. The transition therefore can take place 
extremely fast.) 



7 



That the resonance effect exists should be expected from the physical intuition gained from the classical 
solution. For the special case of an electric field along a fixed direction, say a sinusoidal electric field along 
the direction of xi, the transition probabilities can be calculated by choosing the gauge A(x) = (0, Bx,0) 
and the condition p2 = [H] which in effect makes the geometric operator equal to the identity. For the 
general situation where the electric field changes direction, the nonadiabatic factor can be deduced only by 
first separating out the geometric operator in the time evolution. This includes the cases studied in the 
previous section and in this section. 



8 Appendix A 

To derive the factorization of U(t,0), first switch to the Heisenberg picture. The equations of motion are 

qB ■ qB ■ 

7T = -iu-k + i — R(t), — —— e M „i? M (i), (22) 

where 71 = 71!+ Z7r 2 , and R(t) — R\(t) + iR 2 (t). The solution to the Heisenberg equations can then be 
expressed as 

t 

ir(t) = 7r(0)e-^* + i^-e~ lujt [ e lus ^-R{s)ds, (23) 
c J as 

o 

Wfi (t) = w„(0) - ^e„ u R v {t), (24) 

where we assume R u (0) — 0. The homogeneous terms in the expressions for ir(t) and w^(t) are generated by 
the usual dynamical operator D(t) = exp(—iHot/K). To produce the extra terms in the expression for 7r(t), 
and Wu(t), respectively, using an operator W(t), such that D(t)W(t) recovers the whole solution, it suffices 
for W(t) to satisfy: 

t 

W\t)n(p)W(t) =7r(0)+i— [ e luJs ^R{s)ds, 

c J as 
o 

W\t)w^(0)W(t) = Wfi (0) - ^e^(RAt) - ^(0)). 

In view of the commutation relations (1), which imply [7r,7r'] = 2hqB/c, and from the formula exp(— B) A exp(B) 
A + [A, B] with the condition that [A-B] commutes with A and B, it is clear that W{t) can be chosen to 
be the product of two mutually commuting operators, generated by (1, 7r(0), tt^{0)) and (1, wi(0), W2(0)) 
respectively. Each of theses operators produces a translation for either 7r(0) or uy(0) while leaving the other 
unchanged. Writing W(t) as W(t) = J(t)M(t), we can make the following choice for J(t) and M(t), 

t t 

J(t) = Texp (i^M J e*»±R{8)da + i^- J e~ ius -^j-R* (s)ds\ , (25) 

o o 

M(t) = Pexp(- i ^- 1 i i ; M (0)i? M (i)), (26) 

where Texp stands for time-ordered exponential. It's different from the direct exponential by a numeri- 
cal phase factor only, similar to the path-ordered exponential. Therefore, it can be directly checked that 
D(t)J(t)M(t) recovers the solutions to the Heisenberg equations. 

To verify that D(t)J(t)M(t) not only recovers the solutions to the Heisenberg equations for ir and w^, 
but in fact is the time evolution operator corresponding to H , we now verify that it satisfies the Schrodinger 
equation. Note that M(t) commutes with both D{t) and J(t), so we have 

ih-{D(t)J(t)M(tj) =ih[-(D(t)M{t))]J(t)+ihM(t)D(t)-J(t). 



8 



It is straightforward that 

ih[^(D(t)M(t))]j(t) = (H o (0) + wMt))D(t)J(t)M(t). 

To calculate ihM{t)D{t)§- t J{t), first observe that 

D{t)n{Q)D\t) = Dt(-t)7r(0)D(-i) = 7r(0)e iwi , (27) 
£)(t) 7 r t (0) J D t (i) = (D(t)ir(Q)D\t)y = ^(0)e~ iut . (28) 

Therefore 

d 

ihM(t)D{t) — J{t) = (-n*(0)R(t)/2-%(Q)R*(t)/2)D(t)J(t)M(t), (29) 
= -(7: 1 {0)R 1 {t)+Tr 2 {0)R 2 (t))D{t)J{t)M{t). (30) 
Combining terms and from the definitions of tt^ and w^, we now have 

ih^- t (D(t)J(t)M(t)) = (H + t^-x 1 R 2 -^-x 2 R 1 )D(t)J(t)M(t), (31) 
= H(D(t)J(t)M(t)). (32) 
Therefore, we conclude that the time evolution operator corresponding to H is 

U(t, 0) = D{t)J{t)M{t) = M(t)D(t)J(t). (33) 
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